preprint - UTF 357 



Finite Temperature Effects for Massive Fields 
in D-dimensional Rindler-like Spaces 



Andrei A. BytsenkoQ 

State Technical University, St. Petersburg 195251, Russia 

Guido CognolaQ and Sergio ZerbiniQ 

Dipartimento di Fisica, Universita di Trento 
and Istituto Nazionale di Fisica Nucleare, 
Gruppo Collegato di Trento, Italia 



August 1995 



Abstract: The first quantum corrections to the free energy for massive fields in D -dimensional 
space-times of the form 1R x ]R + x J\A N ^ , where D = N + 1 and M. N ^ is a constant cur- 
vature manifold, is investigated by means of the (^-function regularization. It is suggested 
that the nature of the divergences, which are present in the thermodynamical quantities, 
might be better understood making use of the conformal related optical metric and associ- 
ated techniques. The general form of the horizon divercences of the free energy is obtained 
as a function of free energy densities of fields having negative square masses (absence of the 
gap in the Laplace operator spectrum) on ultrastatic manifolds with hyperbolic spatial sec- 
tion H N ~ 2n and of the Seeley-DeWitt coefficients of the Laplace operator on the manifold 
M. N ~ 1 . Furthermore, recurrence relations are found relating higher and lower dimensions. 
The cases of Rindler space, where M N ~ l = M N ^ and very massive D-dimensional black 
holes, where Ai N ~ l = S N ^ are treated as examples. The renormalization of the internal 
energy is also discussed. 
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1 Introduction 

Recently there has been a renewed interest in the physics of black holes. Several issues like 
the interpretation of the Bekenstein-Hawking classical formula for the black hole entropy, the 
puzzle of loss of information in the black hole evaporation and the interpretation of the Hawking 
temperature have been discussed (see, for example, the review Furthermore in many papers, 
it has been pointed out that it should be desirable to have the usual statistical interpretation 
of the black hole entropy as the number of the gravitational states at the horizon and to try to 
understand the dynamical origin of the black hole entropy (see, for example |^]). 
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On general grounds, the density of levels as a function of the mass M, of a D-dimensional 
black hole should read H (for details, see Sec. 6.3) 

O(M) ~ C D (M) exp (^^M^j . (1.1) 

where Cyj(M) is a quantum prefactor and Gr> is related to the generalized Newton constant 
(see Eq. ( |6.9D ). For a 4-dimensional black hole we have 0] 



Q(M) ~ C 4 (M) exp(4vrGM 2 ) , (1.2) 

4-irGM 2 being the Bekenstein-Hawking classical entropy pi, |(| |7j. As 't Hooft has pointed 
out, the prefactor C^{M) (the first quantum correction to the classical result), which is usually 
computable for quantum fields or extended objects (such as strings or p-branes) in an ultrastatic 
space-time background, turns out to be divergent. This prefactor may be regarded as the 
contribution associated with the first quantum correction to the classical free energy. Several 
different methods have been used in dealing with such divergences, for example "the brick wall 
method" |3], ||, ||], the conical singularity method |l^, 11, 12, [13], |l4|, |15|| , critically discussed in 



1 16, [L7]] and the related "entaglement entropy method" [18, [H], |2( 

The horizon divergences have also been associated with the information loss issue of black 
holes 0, |[ and their physical origin, for quantum fields pll, 22] or strings [23, 24, 25, 26 1, may 



be described by the following simple considerations. 

In a D-dimensional static space-time with horizons, the equivalence principle implies that a 
system in thermal equilibrium has a local Tolman temperature given by T(x) = T j \/—goo( x ), 
T being the asymptotic temperature. Roughly speaking, near the canonical horizon (this means 
that the quantity goo nas simple zeros), a static space-time may be regarded as a Rindler-like 
space-time. We will show that, if one denotes by p the proper distance from the horizon, one gets 
for the Tolman temperature T(p) = T/p. As a consequence, the total entropy for a quantum 
bosonic gas reads (omitting a multiplicative constant) 

//■oo rco 
dx / T(p) - 1 dp ~ AT ' 1 / p- D+1 dp , (1.3) 
JO JO 

where A is the integral on the transverse coordinates x, namely the area of the horizon. The 
latter integral is clearly divergent. Introducing a horizon cutoff parameter e we may rewrite it 
as 

rA A 

S ~ AT p^ 1 dp ~ AT In — , for D = 2 A infrared cutoff 

Je ' £ 

/OO ArpD — X 

p- D+1 dp ~ _ 2 , for ^2. (1.4) 

For the sake of generality, we write down the asymptotic high temperature expansion for the 
entropy of a quantum gas on a D = N + 1-dimensional static space-time defined by the metric 

ds 2 =g 00 ( X )(dx ) 2 + g ij ( X )dx i dxi , x = {^'}, i,j = 1, N , (1.5) 

g(x) denoting its determinant. Again the equivalence principle leads to 

s ^ t n f (jwy N/2 dx N (16) 

As a consequence, the horizon divergences depend on the nature of the poles of the integrand 
(9/900 )~ N ^ 2 - In general, for extremal black holes, where goo has higher order zeroes, the diver- 



gences are much more severe than in the non extremal case (see for example Refs. [27, 17, |2 



These considerations suggest the use of the optical metric g^ u = g^u/doo, conformally related 
to the original one, in order to investigate these issues. It is our opinion that the conformal trans- 
formation techniques are particularly suitable for studying finite temperature effects for fields 
in space-times with horizons and here we would like to present some examples of computation. 



This method has already been appeared for example in Refs. [29, [3(], 31] and has been recently 
used in the horizon divergence problems in Refs. ||22|, ||, [H], 17, [32], 33]. See also [34], where the 
same result is obtained with a different approach. 

One of the purposes of this paper is to implement this idea in the case of massive fields in 
D-dimensional Rindler-like space-times, we are going to introduce. Let us consider static space- 
times admitting canonical horizons and having the topology of the form M x M + x J\A. N ~ l . The 
metric reads 



ds 2 = Tj—dxn + dp 2 + d<r 2 N _ l , (1.7) 

where r# is a dimensional constant, b a constant factor and dafj_ 1 the spatial metric related 
to the N — 1-dimensional manifold M. N ~ X . If M. N ~ l = M 1 ^^ 1 , b = 1 and r# = 1/a, a being 
a constant acceleration, one has to deal with the Rindler space-time. Quantum fields in such 
a case have been considererd in many places, see for example Refs. [B5L [36], E37L 38, 39]. If 



M N ~ l = S N ~ l , b = (D — 3)/2 and is the Schwarzschild radius of a black hole, then we shall 
show that one is dealing with a space-time which approximates, near the horizon and in the 
large mass limit, a D-dimensional black hole (see Sec. |3,2| ). 

It is well known that space-times with canonical horizons admit a distinguished temperature, 
the (Unruh) Hawking temperature. There are several ways to compute it, one of the simplest 
makes use of the relation with the related surface gravity. The other one consists in imposing 
the absence of conical singularities in the Euclidean continuation of the space-time itself (7| . For 



the metric (1.7), one obtains 



0H = ^ ■ (1.8) 

One can arrive at the same result working without using the Euclidean continuation method, 
but making use of the principle of local definiteness in quantum field theory [36, fKJ]. Note that 



these method are no longer equivalent when one is dealing with extremal black holes [41]. It is 
also important to stress that the variable p defined by the metric in Eq. ( |1.7| ) has the meaning 
of radial proper distance between the horizon and a point outside it and so, the divergences of 
thermodynamical quantities are automatically expressed in an invariant way. 

The contents of the paper are the following. In Sec. a review of the necessary conformal 
transformation techniques is presented. In Sec. ||, we consider a Laplace-type operator defined 
on a class of D = N+ 1-dimensional space-times, whose spatial sections have metrics conformally 
related to M N = M + x M N ^. Since in general Ai in a non compact manifold, the Laplace 
operator has a continuum spectrum and a general form for the Plancherel measure, which is 
the analogue of the degeneracy in the case of discrete spectrum, is presented. The measure 
is used in Sec. ^] in order to obtain a useful form for the trace of the heat kernel, which is 
necessary for the derivation of the free energy, which we derive in Sec. ||. It is pointed out that 
in the Rindler case, the spatial section of the conformally related space-time turns out to be 
an iV-dimensional hyperbolic manifold. In this case, the massless scalar field can be treated 
without approximations. In Sec. |6] some applications to the statistical mechanics of a scalar 
field in D-dimensional Rindler and black hole space-times are presented and the divergences of 
the first quantum corrections to free energy and entropy are given. Finally, we end with some 
conclusions in Sec. ^ and with a resume of heat kernel, ^-function and free energy on constant 
curvature manifolds in the Appendix. 



2 Conformal transformation techniques and optical manifold 

In this section we shall briefly summarize the method of conformal transformations using £- 
function regularization |29|, 42, fMf| . These techniques permit to compute all physical quanti- 
ties in an ultrastatic manifold (called the optical manifold and, at the end of calculations, 
transform back them to a static one. This method is particularly useful in dealing with finite 
temperature effects for quantum fields, since these effects can be easily investigated in ultrastatic 
space-times. 

To start with, we consider a non self-interacting scalar field on a D = N + 1-dimensional 
static space-time defined by the metric dOj), i.e. 



ds z = fi( o(x)((ix ) 2 + gij(x)dx l dx J , 



1,...,N. 



The one-loop partition function is given by (we perform the Wick rotation xq 
all the differential operator one is dealing with will be elliptic) 



Z 



exp 



(pL£>(pd D x 



where ^ is a scalar density of wight —1/2 and the operator Ljj has the form 
L D = -A 9 D + m 2 + ^R 9 . 



(2.1) 
-it, thus 

(2.2) 
(2.3) 



Here m (the mass) and £ are arbitrary parameters, while A^, and R 9 are respectively the 
Laplace-Beltrami operator and the scalar curvature of the manifold in the original metric g. 
The ultrastatic metric g^ u can be related to the static one by the conformal transformation 



SV(x) 



„2«t(x) 



(2.4) 



with <r(x) = — ^ In goo- I n this manner, goo = 1 and gij = gij/goo (optical metric). Recalling 
that by a conformal transformation (we remind that cp is a scalar density) 



BP 



-2cr 



one obtains 

L D (j> = 

where (,d = 
have (j)Lr)<j> ~- 



[R 9 



A 9 

D 



2(D - l)A 9 D a — (D - 1)(D - 2)g^d fM ad u a] , 
(D - 2) 2 



(2.5) 



D 



-A° D a 



A 9 D + i D {e 2a R 9 - R 9 ) 



f {-A 9 D + CdRP + e~ 2 °[m 2 + (£ - £ D )R 9 }} 4> , 



(2.6) 



(D — 2)/A{D — 1) is the conformal invariant factor. From the latter equation we 
= 4>Ld4>, where, by definition 



L D = e-°L D e- a = -A 9 D + £ D R 9 + e~ 2a m 2 + (£ - £ D )R 9 



(2.7) 



This means that action I = I by definition. Note that classical conformal invariance requires 
the action to be invariant in form, that is I = I[4>,g], as to say Ld = Lq. As is well known, this 
happens only for conformally coupled massless fields (£ = £d). 
For the one-loop partition function we have 



Z = J[g, g] Z , 



(2.1 



A 



where J[g, g] is the Jacobian of the conformal transformation. Such a Jacobian can be computed 
for any infinitesimal conformal transformation [42]. To this aim it is convenient to introduce a 
family of continuous conformal transformations 



Q q 



e 2 ^g»„ 



gZ = J\detgl„\=e D «°^g 



(2.9) 



in such a way that the metric is g^ u or g^ v according to whether q = or q = 1 respectively. 
Then one gets 



In J[g q ,g q +s q ] = In 



Zq+5q 



Sq 



(Att) d / 2 



J k D (x\L q D )a(x)^d D x , 



(2.10) 



where ki)(x\L q D ), is the Seeley-DeWitt coefficient, which in the case of conformal invariant 
theories, is proportional to the trace anomaly. In general, one has the asymptotic expansion 



TV e~ tLD ~ ^ K n (Li))t 



n-D 
2 



K n (L 



1 



D 



(4tt) d / 2 



k n {x\L D )^d D x. (2.11 



M 



If the manifold is without boundary then K n = for any odd n. The heat kernel coefficients 
are computable and depend on invariant quantities built up with curvature (field strength) and 
their derivatives (see, for example, [f46f). 

The Jacobian for a finite transformation can be obtained from Eq. ( |2.10| ) by an elementary 
integration in q 03]. In particular we have 



In J[g,g\ 



1 



dq J kD{x\L q D )a(x)\f(ft d D x , 



(4vr) D / 2 

and finally, making use of the ^-function regularization, one has 
In Z = In Z - In J[g, g] = ^('(0\L D £ 2 ) - In J\g, g] , 



(2.12) 



(2.13) 



where i is an arbitrary parameter necessary to adjust the dimensions and £' represents the 
derivative with respect to s of the £-function C(s|-^d^ 2 ) related to the operator Lp, which is 
given by Eq. ( ^7|) . 

The same analysis can be easily extended to the finite temperature case [j43| . In fact we recall 
that for a scalar field in thermal equilibrium at finite temperature T = 1/(3 in an ultrastatic 
space-time, the corresponding partition function Zp may be obtained, within the path integral 
approach, simply by Wick rotation r = ix° and imposing a j3 periodicity in r for the field 4>{t, x l ) 
(i = 1, N, N = D — 1) [^7], 48, 49, |50|] . In this way, in the one loop approximation, one has 

exp f- f dT [ 4>L D 0d N x) . (2.14) 



in which 
Id 



0(T,Z i )=</>(T+/3,Z i ) 



8 2 - A N + £ D B? + e~ la m z + {t-t D )Ry = -d 2 + L N . 



(2.15) 



Since the space-time is ultrastatic, by means of the ^-function regularization again one easily 
obtain flBTl 



InZr- 



2 



VVQ(-\\L N ) + (2 - 21n2£)ResC(-^7v) 

d_ p_ 

i'-^o ds y/4nT(s) £ =1 



+ lim 



J2 / t s -V 2 e- n2 P 2 / u Tre-^dt 



(2.16) 



where PP and Res stand for the principal part and the residue of the function. The free energy 
is related to the canonical partition function by means of equation 



ft 



In Z fi 







(InZp -In J[g, g}) =F + F p 



(2.17) 



where Fq represents the vacuum energy, which is given by the first term in Eq. (2.16), while 
F/s represents the temperature dependent part (statistical sum). The entropy and the internal 
energy of the system are given by the usual thermodynamical formulae 



S p = P 2 dpF((3) = fdpFp 



(2.18) 



Up = f3dpF{ft) + F = -ed p In Z p . (2.19) 

In a similar way, one can consider spinor fields. It is sufficient to make use of the following 
formal identity |Q| 

F f {l3) = 2F 2 p-Fp, (2.20) 

where on the r.h.s. the spinor quantities are left understood in the formal expression of the 
bosonic free energy. As a consequence, the horizon divergences of the bosonic sector cannot be 
compensated by the corresponding fermionic ones. 

3 Spectral function for rank 1 Riemannian spaces conformally 
related to M N = M + x M N ~ l 

In many interesting physical cases, the Euclidean optical metric my be written in the form (see 
Eq. flD|) and Sec. 6) 

2 

ds 2 = dT 2 + r -^ [dp 2 + da 2 N _^j , (3.1) 

where dr = bdx°, rjj being a characteristic length (for example the horizon radius), {thI P) 2 the 
conformal factor and da\_-y the metric of a N — 1-dimensional manifold . 

Here we derive the spectral measure of the operator Ljv, as defined by Eq. (|2.15| ), acting 
on scalars in the spatial section of the manifold defined by the metric Eq. (3.1). Using such 



equations (for convenience now we put rjj = 1; in this way all quantities are dimensionless; the 
dimensions will be easily restored at the end of calculations) we easily obtain 

dV = p~ N dpdV N -!, 

L N = -A%-g 2 N + Cp 2 , (3.2) 
A N = p 2 d 2 p -(N-2)pd p + p 2 A N _ 1 , 

where Ajv-i is the Laplace-Beltrami operator on the manifold .M^ -1 and dV^-i its invariant 
measure. We have also set qn = (N—l)/2 and C = m 2 +t;R 9 . It should be noted the appearance 
of an effective "tachionic" mass —Q%, which has important consequences on the structure of the 
^-function related to the operator Ljy. 

In order to study the quantum properties of matter fields defined on this ultrastatic manifold, 
it is sufficient to investigate the kernel of the operator e~ tLN . To this aim, we will search for the 
spectral resolution of the elliptic operator Ljy. Let ^ ra (x) be its eingenfunctions, namely 

L N ^ra(x) = X 2 ^f ra (x) = X 2 f a (x)(j) ra (p) , (3.3) 



where / Q (x) are the (normalized) eingenfunctions of the reduced operator Ljv-i = — Ajv-i + C 
with eigenvalues A 2 . Note that we assume C to be constant. This means that we restrict our- 
selves to consider only manifolds A4 D with constant scalar curvature, or alternatively minimally 
coupled fields. Moreover, to avoid null eigenvalues we suppose C > 0, but the results can be 
easily extended to the case C = 0. Note that the spectrum of Lat_i could also be continuum. 
The differential equation which determines the continuum spectrum turns out to be 



p 2 d p — {N — 2)pd p - p 2 X 2 a + Q N + X r <j> ra (p) = . 



(3.4) 



The only solutions of the latter equation with the correct decay properties at infinity are the 
Bessel functions of imaginary argument Ki r {pX a ) with A 2 , = r 2 (if C = the operator Lat_i has 
a zero mode and gives other solutions to Eq. ( |3.4| )). Thus we have 

(3.5) 



JV-l 

<Pra(,P) = P 2 K ir (p\ a ) . 

If we interpret in the sense of the distribution the following innner product 
dp 



f°° dp f 

*m,*r'a') = J q J dV N ^i^* a (x)^ r > a >(x) , 



P 

we have the normalization condition 

S(r — r') 



(3.6) 



(3.7) 



where p(r) is the spectral measure associated with the continuum spectrum. Thus, for the heat 
kernel of any suitable function /(Ljv) we may write 



POO 

<x\f(L N )\x' >= / drp(r)f(r 2 )J2K a (^ra(x) 
Jo „ 



(3.8) 



The measure p{r) may be determined in the following standard way (Harish-Chandra's method 
which makes use of the asymptotic behaviour of the Mac Donald functions at the origin. 
From Eq. fl3.4| ) and its complex conjugate and making use of Eq. ( |3.5| ) one arrives at 

(A 2 — -^r') 



limp (7V 2) {dp(/)* ra 4>r'a' - KadpCfir'a') 
p^O 



By means of Eqs. ( pTq ) and (3.7) we get 
5{r — r') u 



H{r) 



lim -= 



[d u K* r (u)K ir ,(u) - K* r (u)d u K ir ,(u)} , 



(3.9) 



(3.10) 



where again, the limit has to be understood in the sense of distributions. Recalling that for 
u -» 



Kir (it) 



-IT 



+ T(ir) 



ir g % 



and 



lim 



u^O x 

one finally has 
fi(r) 



^ni5(x) , 



7r|r(ir)| 2 



r sinh irr 



7T^ 



(3.11) 



(3.12) 



(3.13) 



which is in agreement with the 2-dimensional Kontorovich-Lebedeev inversion formula |53|. 



Since our aim is to evaluate the trace of functions of using Eq. (2.11), in particular 
TV exp(— tLjsf), it is convenient to make the sum over a, introducing the total spectral measure 

n N (r,x) =v(r)p N ~ 1 J2\Ux)K tr (p\a)\ 2 (3.14) 

a 

and integrate on the manifold defining 

WW= / Pi N (r,x)dV = p(r) rj2\Kir(pXa)\ 2 ^, (3.15) 
Jm n Jo a P 

where we have used the normalization properties of f a . In this way, for any suitable function 
/(Zjv) we have 

/•CO 

Tr f(L N ) = / f(r 2 )p I (r)dr . (3.16) 
Jo 

As a simple application of Eq. (|3.14| ) let us consider a massless scalar field in a D-dimensional 
Rindler space-time. In this case the optical spatial section turns out to be the hyperbolic space 
H N and the measure pi N should not depend on x, since one is dealing with a homogeneous 
space and it should coincide with the known Plancherel measure. For this case M. N ~ l = M N_1 
and moreover C = 0. The reduced operator Ljv_i = — Ajy_i has a continuum spectrum, the 
eigenvalues being k 2 = k k and the corresponding eigenfunctions /k = (2tt)~^ n ~ 1 ^ 2 exp(ik-x). 
As a consequence 

O N — 1 poo 

<Mr) = /* z » = M (r) j Q \K ir ( P k)\ 2 k N - 2 dk 

\T(ir + 0N )\> 



(47r) JV / 2 r(7V/2) \T{ir)\ 2 



Otv being the volume of the A r -dimensional sphere. Of course, Eq. ( 3.17 ) is the correct Plancherel 
measure of the Laplace operator in H N f54|. 



4 The heat kernel for massive fields 



Here we derive a general expression for /U/(r) by making use of Eq. ( |3.15| ) and then derive the 
trace of the heat kernel, which is needed for the construction of the partition function according 
to Eq. ( 2.16 ). Now we use the Mellin-Barnes representation [55 



Kf r (p\ c 



T{z + ir)T(z - ir)T(z)_ p -2z X -2z dz 



4iv^^Re 2 >i r(z + l/2) 
and observe that, for Rez > (N — l)/2, the sum over a can be done and gives 



(4.1) 



(4.2) 



In integrating over p, one has to pay attention to the fact that the result is formally divergent. 
For this aim we introduce a horizon cutoff parameter e and, when possible, we take the limit 
e — > 0. Then we get 



p(r) 



r(z + ir)T(z-ir)T(z)C(z\L N ^) 



dz . 



(4.3) 



JRez=o(N-i)/2 zT(z + l/2)e 2z 

The integration over z can be done since the meromorphic structure of T-and ^-functions are 
known. In fact, we have [56] 



r(zK(z\L N ^) = £ Kn{L N N _;ll + Js-i(z) , 



(4.4) 



n=0 



D 



where Jn-i is an analytic function. Since the manifold Ai N_1 has no boundary, all K n with 
odd n are vanishing. 

To make the integral we consider the rectangular contour r = { Re z = c, Im z = a, Re z = 
— c, Imz = —a} and observe that the two horizontal paths Imz = ±o give a vanishing contri- 
bution in the limit a — > do, as well as the path Re z = — c in the limit e —* 0. Also the poles in 
the strip — c < Re z < give a vanishing contribution as soon as e — ► 0. Then we have to take 
into consideration only the poles of the integrand in Eq. (|4.3| ) in the half-plane Re z > 0. Such 
a function has simple poles at the points z = 0, z = —n ± zr and z = (N — 1 — n)/2 (n > 0). If 
Z> is even, that is TV" is odd, z = is a double pole. It is clear that all poles with Rez > give 
rise to divergences, the number of them depending on N, while the poles at z = and z = ±ir 
give rise to finite contributions. As a result one obtains 



r jv-2 

I 2 



\- K 2n {L N -i)^N-2n{r) /47T 
W(r) = jV-l-2n b 

| C ; (0|^-i) , C(0|^-i) 



2vr 



2tt 



^?(ir) + ip(-ir) — 2 In ir8(r) 



, (4.5) 



where 



N-2 



is the integer part of the number t ne logarithmic derivative of T and 

6(r) the usual Dirac ^-function. Note that for even N, ^(0|Ljv-i) is vanishing and so the last 
term in the latter equation disappears. 



Now the trace of the heat kernel can be computed by using Eq. ( |3.16 ) with /(r 2 ) = exp(— tr 2 ) 
We write it in the form 



Tr e~ tLN 



N-3 
2 



K 2n (L N -l) K(t\ - A H N-2n ~ Q 2 N _2n) f^_\ ~ 

N - 1 - 2n 

n=0 



+ - 



1 



4Vvrt 



C'(0|L^_ 1 )-2C(0|L w _ 1 )ln- 



2vr 



ijj(ir)e r dr , 



(4.6) 



Tre 



N-2 
2 

E 

n=0 



K 2 „(iw-i)if(t| - A H , 



0N-2n) 



N - 1 - 2n 



4vrV 



+ - 



1 



C(0|L 



AT-lJ , 



(4.7) 



valid for odd and even respectively. Here by K(t\ — A H N-2n — ^ 2 v-2«.) we indicate the 
diagonal heat kernel of a Laplace- like operator on H N ~ 2n . Of course, it does not depend on 
the coordinates since hyperbolic manifolds are homogeneous. Such a kernel is known in any 
dimension [pi]] (see the Appendix). 

As in the previous section, as a simple application of Eq. (4.5), we again consider a massless 

— _ jV— 1 

scalar field in the D-dimensional Rindler space-time. We have Kq(Ln-i) = (47r) ~Vn-i, 
K n = for n > and C( z \^n-i) = for z < (N — l)/2. Here Vjv-i is the volume of the 
manifold M N ~ l (infinite transverse area). Then, using Eq. ( 3.17 ), we immediately obtain 



m(r) = $ N (r)V e 



V, 



N -I 



(4.8) 



which is the integral version of Eq. (3.17). Here V £ may be considered as the volume of H 



N 



A 



5 The thermodynamical quantities 



Now it is quite straightforward to obtain the partition function and then all the others thermo- 
dynamical quantities by means of Eqs. ( gig ) and (|27T7|) . Since the vacuum energy has been 
extensively studied in many papers, here we concentrate our attention on the temperature de- 
pendent part of the free energy (statistical sum) Fp = F@ = — In Zp/f3. Using Eqs. ( (4.6|) and 



(4.7) we get 
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where T^_ 2n indicates the free energy density for a scalar field with (negative) square mass 
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~VN-2n on an ultrastatic manifold with hyperbolic H spatial section, which has been 
studied in detail in Ref. |5l| and is given in the Appendix. 



Some remarks on Eqs. ( |5.1| ) and fl5.2| ) are in order. First of all, it has to be noted that the 
parameter (3 is the inverse of the physical temperature only if b = 1. More generally, before 
to interpret P~ l as the temperature in Eqs. ( |5,lD and (|5.2| ), one has to make the substitution 
P — * b(3. The reason is due to the fact that, in order to write the metric ( |1.7| ) in the form ( |3.1| ). 
we have changed the time coordinate according to r = bxQ. 

Independently on the manifold M N_1 , we see that the (non renormalized) free energy has a 
leading divergence of the kind e~( D ~ 2 ^ proportional to the transverse area Vd-2, since Kq and 



J-pj are always non vanishing. More generally, one has 



D-l 



divergences of the kind e 



~{D-2-2n) 



(depending on the manifold and the operator Lat_i) and, for even D, also a possible logarithmic 
divergence. All these divergences are also present in the expressions of internal energy and 
entropy and their expressions can be obtained by means of Eqs. (|2l9|) and flDlf ). For example 
the internal energy reads 
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where W^_ 2n indicates the free energy density for a scalar field with (negative) square mass 
~&%-2n on an ultrastatic manifold with hyperbolic H N ~ 2n spatial section and Uo(e) is the 
vacuum energy. 
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With regard to the internal energy, we have at disposal a renormalization procedure, which 
is well understood for D = 4. In fact, in Rindler and black hole space-times it is known that the 



renormalized stress-energy tensor is finite at the horizon in the Hartle-Hawking state [57, pi] , 
corresponding to the temperature (3 = (3h- This is equivalent to write 

U{f3) ren = U 6 (p) - U £ (f3 H ) + finite part , (5.5) 

where U £ {j3) is the divergent part of the internal energy and it may be read off the Eqs. ( |5.3|) and 
(5.4). Thus, the divergences are present in the expression of the renormalized internal energy, 



but only for some particular value of /?, say /3h (Ay)- For example, in the case of Rindler space- 
time, such a value is flu = 2-7ra~ 1 , the Unruh temperature (here a is the acceleration), while in 
the 4-dimensional black hole background one has f3jj = 8irMG, the Hawking temperature. 

In the general case, we may use of the same renormalization procedure. Note, however, that 
the corresponding renormalized partition function, free energy and entropy remain divergent 
also at the distinguised temperature f3 = (3h- 



6 Some physical applications 

As simple physical applications of the general formulae derived in Sec. ||, here we consider the 
cases in which M N ~ 1 is a homogeneous manifold with constant scalar curvature k. Of course 
we have the three possibilities M N ~ l = M N ' 1 (k = 0), M N ^ = S N ' 1 (k > 0) and finally 
M N ~ l = H N ~ l (k < 0), but here we only consider in more detail the first two cases. The first 
one corresponds to the conformal treatment of the D-dimensional Rindler space-time, while the 
second appears when one studies the physics of black holes near the horizon. For D = 4, this 
case has been studied in Ref. flBl 



6.1 Statistical mechanics for massive fields in the Rindler D-dimensional 
space-time 

As we have already observed, after a conformal transformation, the spatial section of the Rindler 
space-time is of the kind condidered in the paper. For this special case, the curvature of 
is vanishing {k = 0) and so one easily has 
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where C = m 2 has been put since Rindler is a flat manifold. Now, using Eqs. (5.1) and ( |5.2| ) 
together with the two equations above, we obtain 
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valid for even and odd D-dimension respectively. In Eq. ( |6.2| ) 7 is the Euler constant. The 



functions •7 r ^_ 2n can be computed using Eqs_Q), flPD and flO§D in the Appendix. 



For example, when D = 4, using Eq. ( A.13[) , the result is 
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where the transverse area A = V2 has been introduced to compare the latter formula with well 
known results (see for example Refs. & ETSL BHj). There is agreement in the massless case, 
but not in the massive case, where we also obtain a finite contribution. We conclude this section 
with some remarks on renormalization. As we have seen above, in our formalism, massive scalar 
fields in Rindler space-time can be easily treated, because the optical spatial section turns out 
to be the hyperbolic space H 3 and the harmonic analysis on such a manifold is well known. The 
formulae are particularly simple in the massless case. For example, in 4-dimensions, the total 
free energy may be chosen in the form 
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where /% = 2tt is the Unruh temperature (a = 1). As a consequence, the entropy turns out to 
be 
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and it diverges for every finite /3, but is zero at zero temperature (the Fulling- Rindler state), 
which is correct, since we are dealing with a pure state. At j3 = (3jji corresponding to the 
Minkowski vacuum, we have a divergent entropy proportional to the area, regardless of the 
fact that the Minkowski vacuum is a pure state. This is also to be expected, since an uni- 
formly accelerated observer cannot observe the whole Minkowski space-time. Finally with this 
renormalization prescription, the internal energy should read 
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and this is vanishing and a fortiori finite at (3 = Pjj, as it should be. Furthermore, at = 00, 



namely in the Fulling- Rindler vacuum, it is in agreement with the result obtained in Ref. [Sl|. 



6.2 Statistical mechanics for massive fields in a D-dimensional black hole 
background 

Here we consider in more detail the case in which A^ -1 = S 1 ^^ 1 . To justify this choice from 
a physical view point, first of all we show that, near the horizon, a D-dimensional black hole 
may be approximated by a manifold of this kind and so, the thermodynamics can be derived by 
using the formulae of Sec. [|. 

The static metric describing a D-dimensional Schwarzschild black hole (we assume D > 3) 
reads 

-1 

dr 2 + r 2 dn D „ 2 , (6.8) 
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where we are using polar coordinates, r being the radial one and dQ^-2 the D — 2-dimensional 
spherical unit metric. The horizon radius is given by 
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M being the mass of the black hole and Go the generalized Newton constant. The associ- 
ated Hawking temperature reads (3h = 47rrjj/(.D — 3). The corresponding Bekenstein-Hawking 
entropy may be computed by making use of 
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Thus we have 
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From now on, we put Th = 1- It may be convenient to redefine the radial Schwarzschild 
coordinate r = r{p) by means of the implicit relation 
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and time xq = x' /b, b = (D — 3)/2 in order to have goo = p 2 + 0(p 4 ). In the new set of 
coordinates we have 



ds 2 



1 



„3-D 



(P) 



dx'n + 



1 



„3-D 



(P) 



b 2 1 b 2 p 2 

and finally the optical metric reads 
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From the latter equation we see that, near the horizon p = 0, we can set G(p) = 1 and so 
the optical metric assumes the form considered in previous Sections. In this approximation the 
manifold M N ~ l becomes the unit shpere S 1 ^ -1 . We have 



ds 2 ~ -dx'o 2 + 1 



dp 2 + dO.o-2 



(6.16) 



Such a metric can be considered as an approximation of the one of the black hole in Eq. ( |6.14| ) 
in the sense that, near the horizon, the geodesies are essentially the same for both the metrics. 
The metric ( |6.16| ) can be related to a manifold with curvature BP = — (D — 1)(D — 2) + 0(p 2 ), 
then, according to Eq. (2/7), the relevant operator becomes 
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where now C is a positive constant, which takes into account of mass and curvature contributions 
to this order. Note that since for the original manifold R 9 = 0, £ does not appear in the formulae. 
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This effectively happens if we approximate the metric after the optical transformation has been 
done. More simply, one can put £ = £d in Eq. (|2.7|). 

The discussion for arbitrary D = N + 1 is quite involved even though it may be done, since 
the ("-functions of the Laplace-Beltrami operators on S N ~ 1 are known (see the Appendix). 

As a more explicit example, now we consider a scalar field in a 4-dimensional Schwarzschild 



background. Using these techniques, such a case has been considered in Ref. [58], where we refer 
the interested reader for more details. We have rjj = 2MG, b = 1/2, 
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Then, according to Eq. ( |6.17| ), the relevant operator becomes 
Z 3 = -A3 - 1 + c P 2 , 



(6.18) 



(6.19) 



(6.20) 



where C = m 2 + 1/3 takes into account of the curvature RP = —6 + 2p 2 of the optical manifold. 

Now, directly using Eqs. ( jOg ), (|2.17|) , flOl ) and ( |A.13| ), after the replacement (3 -> /3/2 due 
to the redefinition of the Schwarzschild time (remember that 6 = 1/2), for the total free energy 
we obtain 
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where we have written the Jacobian contribution to the partition function due to the conformal 
transformation in the form Af3j e , and now A 



Airr H is the transverse area of the black hole. 



The (- function related to the operator L 2 on the sphere is given in the Appendix by Eq. ( A.18| ) 
with C = m 2 + 1/3, then C(0|i 2 ) = m 2 . 

The leading divergence, due to the optical volume, is proportional to the horizon area [Q], 
but in contrast with the Rindler case, a logarithmic divergence is also present, similar to the 
one found in Refs. Jl^, This is a feature of even dimensions (see Sec. ||). 

Let us briefly discuss the renormalization of the internal energy in this particular case. 
We recall that one needs a renormalization in order to remove the vacuum divergences. These 
divergenges, as well as the Jacobian conformal factor, do not contribute to the entropy. However 
the situation presented here is complicated by the presence of horizon divergences, controlled 
by the cutoff parameter e. In the 4-dimensional Schwarzschild space-time, it is known that the 
renormalized stress-energy tensor is well defined at the horizon in the Hartle-Hawking state [j57| , 
31], which in our formalism corresponds to the Hawking temperature f3 = 0h- The renormalized 
internal energy reads (the dots stay for finite contributions at the horizon, which we do not write 
down because their value depend on the approximation made) 
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which has no divergences for P 
8tt 2 A Alne 



Ph, the Hawking temperature, while the entropy 
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also for P = Ph contains the well known divergent term proportional to the horizon area 
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and, according to Ref. 12], a logarithmic divergence too. Eq. (4.4) is vanishing in the Boulware 
vacuum corresponding to f3 = oo. 
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From this renormalization procedure we get for the renormalized black hole free energy 
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In the general case, the discussion is quite similar and it can be performed by using Eqs. ( |5.1| ) 
or (O) with the replacement (3 ^ {D - 3)/J/2. 



7 Conclusions 

In this paper the first quantum corrections to the thermodynamic quantities of fields in a D- 
dimensional Rindler-like space have been investigated making use of conformal transformation 
techniques and (^-function regularization. In this way, we have worked within the so called optical 
manifold, which is ultrastatic, and the use of finite temperature methods is quite straightforward. 

The general form of the horizon divercences of the free energy has been obtained as a function 
of free energy densities of fields having negative square masses (absence of the gap in the Laplace 
operator spectrum) on ultrastatic manifolds with hyperbolic spatial section H N ~ 2n and of the 
Seeley-DeWitt coefficients i^n^TV-i) of the Laplace operator on M. N ^. Since there exists 
recurrence relations for free energy densities (see the Appendix), it is sufficient to study the 
cases D = 3 and D = 4 [D = 4 and D = 5 for applications to black holes). The leading 
divergence can be seen to be given by the volume of the spatial section of the optical manifold. 
A finite contribution is also obtained and this depends on £(0|Ltv-i) and on its first derivative. 
For D = 4, our results are consistent with the ones obtained with brick wall, path-integral and 



canonical methods [0, 17 



With regard to physical applications, we have used the general results on finite temperature 
field theory in order to investigate the quantum corrections to the Bekenstein-Hawking entropy 
for massive fields in a large mass black hole background. This approach gives rise to a leading 
divergence for the entropy similar to the one obtained for the Rindler case background, but in 
this case other divergent contributions are present and their structure depend on the dimension 
of the space-time considered. Here we have shown how it is possible to get the general form 
valid for an arbitrary dimension and we have explicitly considered the case D = 4. 



We also would like to mention the results obtained in Ref. [63], where the contributions to 
the 4-dimensional black hole entropy due to modes located inside and near the horizon have 
been evaluated using a new invariant statistical mechanical definition for the black hole entropy. 
The finite contributions, namely the ones indipendent on the horizon cutoff, are compatible with 
our results. 

As far as the horizon divergences are concerned we recall that they may be interpreted 
physically in terms of the infinite gravitational redshift existing between the spatial infinity, 
where one measures the generic equilibrium temperature and the horizon, which is classically 
unaccessible for the Schwarzschild external observer. Furthermore, we have argued that they are 
absent in the internal energy at the Unrhu-Hawking temperature. However, they remain in the 
entropy and in the other thermodynamical quantities, as soon as one assumes the validity of the 
usual thermodynamical relations. For D = 4, a possible way to deal with such divergences has 



been suggested in Refs. [Q, 63, 64], where it has been argued that the quantum fluctuations at 
the horizon might provide a natural cutoff. In particular, choosing the horizon cutoff parameter 
of the order of the Planck length (e 2 ~ G), the leading "divergence", evaluated at the Hawking 



temperature, turns out to be of the form of the the "classical" Bekenstein-Hawking entropy. 
This seems a reasonable assumption, because we have worked within the fixed background 
approximation. However one should remark that other terms are present, giving contributions 
which violate the area law. A more elaborate discussion for D = 4 can be found in Ref. fl6E}| . 
Alternatively, one may try to relate the horizon divergences to the ultraviolet divergences of 
quantum gravity, thus arriving at the theory of superstring progagating in a curved space-time 
H or at the renormalization group approach [p6| . 

Finally, we mention that there has been the proposal to absorbe the horizon divergences, at 
least for D = 2, 3, 4, by making use of the standard ultraviolet gravitational constant renormal- 
ization H |67|, 68, 69, 70 1 . This proposal is essentially based on the use of Euclidean section with 



a conical singularity and associated heat kernel expansion. The problematic issue consisting 
in dealing with a finite temperature theory in a non ultrastatic space-time is solved working 
within a non vanishing conical singularity and interpreting the deficit angle of the Euclidean 
compactified time as the inverse of the equilibrium temperature (the absence of the conical 
singularity gives the correct Hawking temperature). However, the resulting partition function 
has, apparently, a wrong dependence on this "temperature". Furthermore, it seems that the 
only divergences present are the usual ultraviolet ones associated with the definition of the par- 
tition function as determinant of an elliptic operator. These divergences are then absorbed in 
the gravitational constant renormalization. However, the naive use of (-function regularization 
should get rid off these ultraviolet divergences. Thus it seems to exist a disagreement between 
this approach and our approach based on the conformal transformation techniques. It is our 
opinion that this disagreement might depend on a non commutative property present in the 
evaluation of heat kernel trace on a cone. It should be interesting to elucidate this issue. 
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A Heat kernel, (-function and free energy in constant curvature 
spaces 

To start with, let us write down a representatation for the statistical sum which may be useful 
for investigating the high temperature expansion. For detailed derivations we refer the reader 



to Ref. [51], where all material of this Appendix can be found. It reads 



T K = ~- L dz C( z )R^ ~ l)((^\A N )(3- z , (A.l) 
TnjRez=c 2 

where An = — Ajv + a 2 + kq\, acting on fields in H N or SV, with a a constant which may go 
to zero. 

A.l Hyperbolic manifolds 

Here we only report the equations concerning heat kernel, ("-function and free energy for scalars 
fields in hyperbolic manifolds, we need in the paper. 

The measure for the Laplace-Beltrami operator acting on scalar fields in H N satisfy the 
recurrence relations 
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which permit to derive recurrence formulae for all others quantities one is interested in. In 
particular, for the operator An (here k = — 1) we have 
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where K(t\AN) and ((s\An) are densities. Finally, for the free energy density, one gets 
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Using equations above, one obtains the desired quantities for fields in H , starting from H 3 or 
H 2 according to whether N is odd or even. For the space we have 
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In Eq. ( [AlOD , B n are the Bernoulli numbers and the series is convergent for < t < 2ir. 

Note that in the paper the relevant operator concerning hyperbolic metric is An_ 2ti with 
a = 0. In this limit, equations above notably simplify, especially in the odd N cases. For 
example, from Eq. ( |A.8| )in the limit a^Owe have 
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A. 2 Spheres 



As has been shown in Ref. 51], the recurrence relations ( [A.3j - A~5| ) are also valid for manifolds 



with positive constant curvature (k = 1). This means that for spheres, heat-kernel, £- function 
and free energy can be obtained by the knowledge of the same quantities for S 1 and S 2 . In 
the paper we only need heat-kernel coefficients and ^-function. F or the circle S 1 we have the 
representation 
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where T is an open path in the complex plane going (clockwise) from oo to oo around the 
positive real axis enclosing the point z = Q2 = 1/2 (see Ref. |5l| ). Note that in the paper the 
operator concerning spheres is L/y-i, so we have to choose a 2 = C — £?jv-i- ^ course, other 



representations are at disposal and can be found for example in Refs. [[Q], 54, 51]. For example 
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whence ((0\L 2 ) = C - 1/3. 
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